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(57) Abstract: The invention relates to a method of de- 
termining the relative position of a mobile element in re- 
lation to the known position of a reference station. The 
method of calculating the relative position comprises the 
following steps, namely consisting in: a) selecting a lin- 
ear combination of satellite transmission frequencies LI 
and L2 from a predetermined list comprising at least two 
linear frequency combinations; b) calculating a precise 
relative position Pp of the mobile in relation to the refer- 
ence station from the linear combinations of pseudo-dis- 
tance corresponding to said linear combination and an 
estimated position Pe of the mobile element in relation 
to the reference station; c) selecting the next linear com- 
bination from the list if it exists and, in this case, repeat- 
ing step b) using the precise position Pp as the estimated 
position; and d) repeating step c) for all of the linear 
combinations from the list. 

(57) Abrege : L* invention concerne un proc&ie' de 
determination de la position relative d'un mobile par 
rapport a la position connue d'une station de reference. 
Le calcul de position relative comprend, les Stapes 
consistant a - a) choisir une combinaison lineaire des 
frequences LI et L2 d' Emission de satellites parmi 
une liste preVJ^terminee comportant au moins deux 
combinaisons lineaires de frequences,- b) calculer 
une position relative precise Pp du mobile par rapport 
a la station de reference, a partir des combinaisons 
lineaires de pseudo-distances correspondant a ladite 
combinaison lineaire et d'une position estimee Pe du 
mobile par rapport a la station de 

[Suite sur la page suivante] 
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Publiee : 

— avec rapport de recherche Internationale 

— avant I 'expiration du delai prevu pour la modification des 
revendications, sera republiee si des modifications sont 
recites 



En ce qui concerne les codes a deux lettres et autres abrevia- 
tions, se referer aux "Notes explicatives relatives aux codes et 
abreviations" figurant au debut de chaque numero ordinaire de 
la Gazette du PCT. 



reference,- c) choisir dans la liste la combinaison lingaire suivante si elle existe, et r&terer dans ce cas Tetape b) en considerant 
comme position estimee, ladite position precise Pp, - d) reiterer Tetape c) pour toutes les combinaisons lineaires de la liste. 
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PROCEDE ET DISPOSITIF DE DETERMINATION DE LA POSITION 
RELATIVE DE DEUX POINTS 



L'invention concerne la determination precise de la position 
relative de deux points qui peuvent etre 6loignes de plusieurs dizaines de 
5 kilometres, a partir de signaux de positionnement par satellites. 

Le domaine d'application est celui des techniques necessitant de 
connattre en trois dimensions et avec une precision centimetrique, la position 
d'un mobile par rapport a la position connue d'une station de reference 
eloignee de plusieurs dizaines de kilometres. On peut citer la geodesie, la 

10 topographie, I'hydrographie 

Pour la determination de la position relative d'un mobile par 
rapport a une station de reference, on utilise couramment des moyens de 
mesure de position par satellites, utilisant par exemple les signaux radio emis 
par les satellites du systeme GPS (Global Positioning System) ou d'autres 
1 5 systemes analogues (systeme GLONASS, futur systeme GALILEO). 

Dans le systeme GPS, le signal emis par un satellite est code et 
on utilise le temps mis par ce signal pour atteindre le point a localiser pour 
determiner la distance entre ce satellite et ce point, de preference appelee 
pseudo-distance pour tenir compte d'erreurs de synchronisation entre 
20 I'horloge du satellite et celle de la station. Ces erreurs de synchronisation 
sont classiquement 6liminees par calcul des lors qu'on recoit les signaux 
d'au moins quatre satellites differents. La determination de la distance entre 
le point a localiser et plusieurs satellites permet, connaissant les 
coordonnees g6ographiques des satellites, de calculer les coordonnees du 
25 point a localiser, le plus souvent des coordonnees exprimees en latitude, 
longitude et altitude dans un repere terrestre fixe. 

Pour determiner la position relative d'un mobile par rapport a une 
station de reference, on utilise une methode dite « GPS differentiel » qui 
consiste a localiser un point par rapport a une station de reference et non par 
30 rapport a un repere terrestre independant : en mettant en place un recepteur 
a la station de reference, on peut determiner, en utilisant les mesures faites a 
la station et au mobile, la position relative du mobile par rapport a la station 
de reference. 
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L'avantage de cette methode est qu'elle permet d'accroftre la 
precision du positionnement. En effet les distorsions des mesures liees aux 
aleas de la propagation des signaux radiosatellitaires sont le plus souvent 
fortement correlees dans I'espace, et la connaissance precise de la position 
5 de la station de reference permet, en comparant les mesures faites a la 
. station aux distances theoriques, de les compenser en grande partie. 

Le temps de propagation est determine d'une part en reference a 
un instant de repere du code pseudo-aleatoire qui module une frequence 

10 porteuse emise par le satellite, cet instant de repere du code permettant 
notamment de determiner la position approximative du mobile c'est-a-dire 
avec une precision de quelques metres a quelques dizaines de metres ; le 
temps de propagation est determine d'autre part en reference a la phase de 
la porteuse recue, la mesure sur la phase, moins bruitee que la mesure sur 

15 le code, permettant de determiner plus precisement la position du mobile 
C ' es t-a-dire avec une precision centimetrique a la condition toutefois qu'on 
leve rambigu'fte sur le nombre de tours de phase puisque la phase ne peut 
etre a priori connue qu'a 2n pres, 2n correspondant a une distance egale a la 
longueur d'onde du signal radiofrequence emis par les satellites. 

20 Dans ce qui suit on ne s'interessera qu'aux mesures de phase, les 

mesures de position de code pouvant etre faites d'une maniere classique. On 
considerera done que les pseudo-distances fournies par le recepteur GPS du 
mobile ou de la station de reference sont essentiellement des valeurs 
numeriques de phase, une valeur de phase etant directement convertie en 

25 valeur de distance, connaissant la longueur d'onde du signal radio emis par 
les satellites. 

Le point central des techniques de positionnement centimetrique 
utilisant les mesures de phase est le calcul preliminaire dit « de 
I'initialisation » dans laquelle on resout le probleme des ambiguTtes sur le 

30 nombre de longueurs d'onde. Ce calcul necessite classiquement la 
connaissance prealable d'une position estimee du mobile, qui peut etre 
obtenue en particulier par une methode telle que celle decrite dans les 
brevets FR 2 715 230 et FR 2 764 708. Cette position estimee est ensuite 
recalee sur la position precise, puis validee au cours de ce calcul 

35 d'initialisation. 
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On considdre plus particulidrement dans la suite P6tape de 
recalage de la position estimee vers une position precise. 

La qualite de cette position precise depend notamment de 
l'6loignement du mobile par rapport & la station de reference. 

5 En effet, les imperfections de la nfethode differentielle viennent 

d'abord du fait que les signaux radio-satellitaires ne rencontrent pas 
exactement les memes conditions de propagation sur les trajets satellite- 
station et satellite-mobile. Les differences des conditions rencontres, a peu 
pr&s nulles k proximite immediate de la station, s'accroissent naturellement 

1 o avec P6loignement. 

Cette difference est principalement due a Pionosphere qui est 
travers6e par les signaux satellite-station et satellite-mobile en des points 
differents alors que Pionosphere n'est pas un milieu homogene. Les mesures 
differentielles bas6es sur les temps de propagation des signaux satellite- 

15 station et satellite-mobile sont alors affectees par cette difference. Cette 
difference peut conduire a une erreur de position du mobile par rapport k la 
station de reference pouvant atteindre 1 k plusieurs cm par km 
d'eloignement. Ainsi, pour un 6loignement entre la station et le mobile 
superieur a une distance de Pordre de 1 0 km, on ne peut garantir la position 

20 du mobile par rapport a la station de reference avec une precision 
centinrfetrique. 

Une premiere solution d6crite dans le brevet n° 2 764 708 propose 
d'une part de r£duire le temps de calcul d'initialisation en particulier le temps 
25 de calcul d'une position inambigue approchee en utilisant notamment des 
combinaisons lineaires des frequences demission L1 et L2 des satellites du 
sysfeme GPS. Elle propose d'autre part de reduire Perreur ionosph6rique ; la 
reduction de Perreur ionospherique s'applique au cours de Petape de 
recalage. Elle consiste a calculer k partir de la position inambigue 
30 approchee, d'une part une position (XL1, YL1, ZL1) pour L1 et d'autre part 
une position (XL2, YL2, ZL2) pour L2, la position precise (X, Y, Z) resultant 
alors de la combinaison lin^aire suivante : 
X = (1,65 XL1-XL2)/0,65 
Y = (1,65 YL1-YL2)/0,65 
35 Z = (1 ,65 ZL1 -ZL2)/0,65. 
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Mais du fait du calcul d'une position sur L1, entache d'une erreur 
ionospherique E et du calcul d'une position sur L2 entache d'une erreur 
ionospherique 1 ,65*E, il n'est pas toujours possible de lever les ambiguites 
lorsque I'erreur ionospherique s'accroit ce qui se produit lorsque la distance 

5 entre le mobile et la station de reference s'accroit. 

Une autre solution classiquement propos6e consiste & mettre en 
place non plus une, mais plusieurs stations de reference constituant ce qui 
est communement appeie un « Roseau ». Selon cette technique, il est 
possible de connaTtre, non seulement les erreurs mesurees en un point 

10 comme on le fait dans le cadre du « GPS differentiel », mais aussi leur 
gradient devolution dans la zone. On compense done en grande partie Feffet 
des d£corellations spatiales des erreurs. Cette solution est efficace, mais elle 
est bien sOr lourde et coOteuse k mettre en ceuvre en raison de 
Infrastructure qu'elle necessite et du coOt des communications entre 

15 stations et mobile. De plus, une telle infrastructure n'existera pas partout. 

II a egalement 6te propose une methode bas6e sur I'exploitation 
du fait que I'erreur ionospherique est fonction de la frequence (en 1/f 2 en 
premiere approximation). 

II est alors possible de determiner cette erreur ou de la r6duire, 

20 voire de I'eiiminer, en remplagant dans les calculs la frequence f (designee 
par L1 ou L2 dans le cas du systdme GPS) par une combinaison Iin6aire des 
frequences porteuses des signaux 6mis par les satellites, c'est-&-dire par une 
combinaison lin§aire de L1 et L2. 

Le resultat d'une combinaison lineaire de L1 , L2 est une nouvelle 

25 frequence L3 a laquelle correspond une longueur d'onde dite longueur 
d'onde apparente. 

Par exemple, dans le cas du systeme GPS pour lequel L1 = 
1,57542 GHz (ce qui correspond k une longueur d'onde d'environ 19 cm) et 
L2= 1,22760 GHz (ce qui correspond a une longueur d'onde d'environ 24 

30 cm), la combinaison de frequences dite « lono Free », 9L1-7L2 permet 
d'eiiminer presque comptetement I'erreur ionospherique. La longueur d'onde 
apparente correspondante est de 5 cm. 

Cependant cette methode est tr&s difficile k appliquer du fait de la 
grande difficulte k lever les ambiguTtes sur des longueurs d'onde aussi 

35 courtes. 
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Le but de I'invention est done de proposer un proc6d6 et un 
dispositlf permettant d'obtenir avec une precision centimetrique la position 
d'un mobile par rapport a une station de reference eventuellement 6loign6e 

5 de plusieurs dizaines de kilometres. 

Pour atteindre ce but, I'invention propose un proced6 de 
determination de la position relative d'un mobile par rapport k la position 
connue d'une station de reference, utilisant chacun une antenne de reception 
de signaux radio issus d'une constellation de satellites de positionnement 

10 emettant sur au moins deux frequences L1 et L2, ce proced6 comprenant la 
determination p6riodique pour chacune desdites frequences, d'un jeu de 2p 
pseudo-distances, k savoir p pseudo-distances entre le mobile et les p 
satellites et p pseudo-distances entre la station de reference et les p 
satellites, la fourniture des pseudo-distances a un organe de calcul de 

15 position, et le calcul par cet organe d'une position relative du mobile par 
rapport k la station de reference k partir d'une part des pseudo-distances et 
d'autre part d'une position estim6e Pe du mobile par rapport k la station de 
reference, ce procede etant principalement caracterise en ce que le calcul de 
position relative comprend, pour un jeu donne de 4p pseudo-distances regu 

20 par I'organe de calcul les etapes suivantes consistant a : 

- a) choisir une combinaison lineaire aL1+bL2 desdites frequences 
L1 et L2 parmi une liste pr6determinee comportant au moins deux 
combinaisons lineaires de frequences, 

- b) calculer les combinaisons lineaires de pseudo-distances 
25 correspondant k ladite combinaison lineaire et, a partir de ces combinaisons 

lineaires de pseudo-distances et de la position estimee Pe, calculer une 
position relative precise Pp du mobile par rapport k la station de reference, 

- c) choisir dans la liste la combinaison lineaire suivante si elle 
existe, et reit6rer dans ce cas retape b) en consid6rant comme position 

30 estimee, ladite position precise Pp, et en utilisant le meme jeu de 4p pseudo- 
distances, pour obtenir une position relative encore plus precise, 

- d) r6iterer l'6tape c) pour toutes les combinaisons lineaires de la 

liste. 

On execute ainsi successivement, k partir du meme jeu de 
35 mesures de pseudo-distances plusieurs calculs de position du mobile en 
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utilisant differentes combinaisons lineaires de frequences, la position estimee 
en entr6e d'un calcul etant la position calcuiee & l'6tape precedente. 

Dans la technique anterieure un seul calcul etait effectue et la 
combinaison lingaire utilisee dans ce calcul portait sur des positions 

5 respectivement calcuiees pour chacune des frequences L1 et L2. 

Une caracteristique importante de I'invention est que les 
combinaisons lineaires de la liste sont d£termin£es de telle sorte que les 
longueurs d'onde correspondantes diminuent progressivement et que la 
sensibility aux erreurs ionospheriques diminue egalement progressivement 

10 et de fa$on plus rapide que la longueur d'onde. 

En d'autres mots, la premiere combinaison lineaire de frequences 
est choisie pour que sa longueur d'onde soit grande, afin de faciliter le lever 
d'ambiguTte sur le nombre de tours de phase ; mais en contrepartie cette 
premiere combinaison lineaire peut correspondre a une sensibility importante 

15 aux erreurs ionospheriques. La position estim6e Pe ayant ete am6lior6e suite 
k ce premier recalage, la deuxieme combinaison lineaire de la liste 
correspond k une longueur d'onde plus faible et a une erreur ionospherique 
encore plus faible. Et ainsi de suite, en utilisant des longueurs d'onde et des 
sensibilit6s aux erreurs ionospheriques de plus en plus faibles. 

20 Dans le cadre du systeme GPS, la premiere combinaison peut 

etre la combinaison L1-L2 (a=1, b=-1) ; la demiere peut etre 9L1-7L2 (a=9, 
b=-7), combinaison connue pour etre pratiquement insensible aux erreurs 
ionosph6riques. 

Les combinaisons intermediaires sont de preference les suivantes 
25 (dans I'ordre) : 

2L1-L2 ; 3L1-2L2 ; 4L2-3L1. 

Selon un autre aspect caracteristique de I'invention, l'6tape b) du 
calcul est effectuee soit en une seule etape utilisant directement les p 
satellites, soit en deux etapes dont la premiere utilise seulement un nombre 
30 reduit p' (p'<p) de satellites et la deuxieme utilise les p satellites. De 
preference, c'est seulement lors de I'utilisation de la premiere combinaison 
lineaire (de longueur d'onde la plus grande) que le calcul est fait en deux 
etapes, les autres combinaisons lineaires etant utilisees selon un calcul en 
une seule etape avec les p satellites. 
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Lorsque Tetape b est effective en deux etapes, elle comprend 
avantageusement les etapes suivantes consistant & : 

b1) calculer une position relative approchee Pa du mobile par 
rapport k la station de reference, a partir de la combinaison lineaire choisie, 

5 de Pe et d'un sous-jeu de 4p' pseudo-distances correspondant a p 1 satellites 
0C1 p' est inferieur a p et 0C1 les p' satellites choisis dans la constellation de p 
satellites sont ceux qui ont, compte-tenu de la ggometrie actuelle de la 
constellation, la plus faible sensibilite a une erreur de la position estimee, 

b2) calculer une position relative precise Pp du mobile par rapport 

10 a la station de reference & partir de ladite combinaison lineaire, de Pa et du 
jeu complet de 4p pseudo-distances. 

Selon une autre caracteristique de I'invention, les 2p pseudo- 
distances entre les satellites et la station de reference sont determines par 
la station de reference et envoyees par radio au mobile qui comporte alors 

15 des moyens de reception pour recevoir ces pseudo-distances et des 
informations de datation de la mesure de ces pseudo-distances. 

Enfin, Invention a pour objet non seulement le proc6d6 de 
determination de position relative d'un mobile par rapport a la position 
connue d'une station de reference dont les grandes lignes viennent d'etre 

20 d£crites, mais aussi un dispositif de determination de la position d'un mobile 
par rapport a une station de reference, capable de mettre en ceuvre ce 
precede . Le dispositif selon invention comprend au moins, dans le mobile, 
des moyens pour recevoir des signaux de positionnement par satellite, et des 
moyens pour recevoir un jeu de 2p pseudo-distances emises par la station 

25 de reference et representant les pseudo-distances entre la station de 
reference et p satellites pour au moins deux frequences porteuses differentes 
L1 et 12, des moyens de determination periodique d'un jeu de 2p pseudo- 
distances entre le mobile et les p satellites, des moyens de fourniture des 4p 
pseudo-distances & un organe de calcul de position, des moyens de 

30 stockage d'une liste de combinaisons lineaires des frequences des porteuses 
des signaux de positionnement, des moyens pour effectuer, h partir d'un 
meme jeu de 4p pseudo-distances, des calculs successifs de position 
relative du mobile par rapport k la position de la station de reference, chaque 
fois a partir d'une combinaison lineaire diff6rente de frequences choisie dans 

35 la liste, d'une position estimee Pe et du jeu de 4p pseudo-distances, la 
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position estirrfee lors d'un calcul avec une combinaison lineaire donnee de la 
liste etant la position relative. calcufee k partir de la combinaison Iin6aire 
prec6dente de la liste. 

5 D'autres caracferistiques et avantages de I'invention apparaTtront 

a la lecture de la description detailfee qui suit, faite k titre d'exemple non 
limitatif et en reference aux dessins annexes dans lesquels : 

- la figure 1 represente sch6matiquement une repartition des 
positions de la station de reference, du mobile et des positions estimee, 

10 approchee et precise du mobile par rapport a la position de la station de 
reference ; 

- la figure 2 represente schematiquement le dispositif selon 
I'invention ; 

- la figure 3 represente un organigramme des calculs effectues. 

15 

On a represents figure 1, la position R de la station de reference, 
eioign6e d'une distance D qui peut etre de plusieurs dizaines de kilometres 
d'un mobile dont la position vraie est M. Sont egalement representees sur 
cette figure les positions estirrfees Pe, approchee Pa et precise Pp du mobile 

20 auxquelles il sera fait reference plus loin. Le point Pe correspond k la 
position M du mobile estimee k quelques metres pres ( 1 a 2 metres par 
exemple), le point Pa & celle du mobile calcufee avec une erreur de un a 
quelques decimetres, le point Pp k celle du mobile calcufee avec une erreur 
maximum de quelques centimetres. 

25 La station de reference et le mobile sont respectivement equip6s 

d'une antenne de reception 10 et 12 de signaux radio issus des satellites de 
positionnement (systeme GPS ou autre) et de moyens de demodulation et de 
traitement des signaux regus. Le mobile calcule periodiquement des pseudo- 
distances entre sa position et la position des satellites k un instant donne. La 

30 station calcule de la meme maniere des pseudo-distances (mesur6es au 
meme instant, ou ramen6es au meme instant de mesure) entre sa position et 
la position des satellites. Comme on I'a dej& vu dans le preambule, les 
pseudo-distances sont classiquement foumies sous forme d'une premiere 
valeur numerique qui d6finit la position temporelle du code pseudo-afeatoire 

35 6mis par un satellite k un moment donne, et d'une deuxieme valeur 
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numerique qui d6f init la phase du signal radiof r^quence module par ce code 
au meme moment. 

Dans les r§cepteurs (mobiles ou stations de reference) simples, k 
faible precision, un seul signal radiofrequence code est utilise par les circuits 

5 de reception. Dans les r6cepteurs plus precis, auxquels la pr6sente invention 
s'interesse, au moins deux signaux de frequence porteuse differente, L1 et 
L2, en provenance des satellites, sont trait6s par le r6cepteur. Le r6cepteur 
calcule done periodiquement des pseudo-distances pour chacune des 
frequences L1 et L2. 

10 II en r£sulte que pour une mesure de position k un instant donn6, 

a partir de p satellites, le calcul de position va faire intervener un jeu de 4p 
pseudo-distances, done 4p mesures de phase, qui sont respectivement : 

- p mesures de phase entre la station de reference et p 
satellites, k la frequence L1 , 

15 - p mesures de phase entre le mobile et les p satellites, a la 

frequence L1 encore, 

- p mesures de phase entre la station de reference et les p 
satellites k la frequence L2, 

- et p mesures de phase entre le mobile et les p satellites, a la 
20 frequence L2. 

La station de reference est 6quip6e de moyens pour 6mettre vers 
le mobile les valeurs nurneriques qu'elle a dSterminees. Le mobile est equips 
de moyens pour les recevoir (antenne 11 sur la figure 2), de sorte que pour 
un calcul de position & un moment donne, le mobile dispose non seulement 
25 de ses propres mesures de pseudo-distances mais aussi des mesures de 
pseudo-distances d6termin§es par la station de reference. 

Dans le mobile, les signaux provenant de la station de r6f6rence 
et les signaux directement regus par le mobile sont utilises par un ensemble 
6lectronique 14 que comporte le mobile et represents figure 2. On pourrait 
30 cependant imaginer la solution inverse oD les calculs sont faits par la station 
de reference, le mobile envoyant k celle-ci les pseudo-distances qu'il a 
calcul6es. Dans ce mode, on connaTt k chaque instant k la station Tendroit 
ou se trouve le mobile, ce qui peut etre interessant dans certaines 
applications. 
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Cet ensemble electronique 14 comprend tout d'abord un circuit de 
reception de positionnement par satellite 16 capable de determiner des 
mesures de pseudo-distances entre son antenne et plusieurs satellites. 

Le circuit de reception 16 fournit p6riodiquement (par exemple 
5 toutes les 100 millisecondes) un jeu de 2p pseudo-distances lorsque p 
satellites sont en vue directe avec Tantenne 12. I! les transmet k un organe 
de calcul de position relative, 18, qui regoit par ailleurs le jeu de 2p pseudo- 
distances correspondant aux 2p pseudo-distances regues des memo p 
satellites par Tantenne 10 de la station de reference et transmises au mobile 
10 par Tantenne 11. 

L'organe de calcul de position relative 18 a pour fonction de 
determiner la position precise Pp de Tantenne 12 du mobile par rapport k 
Tantenne 10 de la station de reference. 

Uorgane de calcul 18 est programme pour effectuer les calculs 
15 desires et est relie aux organes peripheriques necessaires selon les 
applications : afficheur 20, clavier 22, lecteur/enregistreur de donnees ou de 
programmes 24, moyen de transmission filaire ou radio vers un utilisateur, ou 
simple interface de sortie vers une ligne de transmission du resultat des 
calculs. 

20 Les moyens de calcul pr6vus dans Torgane de calcul sont aptes k 

effectuer les operations qu'on va detainer ci-aprds, k partir d'un jeu de 4p 
pseudo-distances regues k un instant donn6 par le mobile. 

Le calcul de position est differentiel, c'est-a-dire qu'on determine 
la position du mobile par rapport a la station de reference (que le calcul soit 

25 effectue par le mobile ou par la station). On peut done considerer dans la 
suite que le calcul de position consiste a calculer la position du mobile a 
partir de mesures differentielles de pseudo-distances, en faisant des 
differences entre pseudo-distances mesurees k la station de reference et au 
mobile. Le calcul fait intervenir une notion de doubles differences de pseudo- 

30 distances (differences entre couples de satellites) que Ton precisera plus 
loin. 

Globalement, en utilisant la notion classique de double difference, 
le principe de calcul differentiel est le suivant : 

- on calcule les positions des satellites k Tinstant de mesure t, 
35 grace aux 6ph6merides des satellites ; 
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- on determine, k partir des mesures differentielles de pseudo- 
distances, les distances differentielles entre la station de reference et le 
mobile selon les axes de vis6e des satellites. On obtient des distances qu'on 
peut consid6rer globalement comme les projections, Ie long de ces axes, de 

5 la distance D entre le mobile et la station de reference ; ce sont des 
distances mesurges ; 

- on calcule pour le m§me instant de mesure les distances le long 
des memes axes entre la station de reference et une position estim6e Pe du 
mobile ; ce sont des distances estinrfees ; 

10 - on determine selon chaque axe la difference entre la distance 

mesur6e et la distance estim£e, qu'on peut appeler grandeur d'£cart ou 
« innovation » selon cet axe ; 

- on calcule k partir de ces grandeurs d'ecart, k I'aide de la matrice 
des cosinus directeurs representant les directions des axes de vis6e des 

15 satellites, les hearts entre la position mesur§e du mobile et la position 
estim§e ; 

- on rajoute k la position estinrfee initiate les hearts calculus, et on 
obtient une position calculee du mobile qui est soit une position definitive Pp, 
soit une nouvelle position estim6e Pe en vue d'une etape ulterieure de calcul 

20 comme on le verra plus loin. 

Comme on Pexpliquera plus loin, les doubles differences 
consid6r6es dans la presente invention sont elabor6es, non pas k partir des 
frequences directement mesurees, mais a partir de combinaisons lin£aires 
des pseudo-distances mesurees a la frequence L1 et des pseudo-distances 

25 mesurees a la frequence L2. 

Un Element essentiel de la presente invention feside dans le fait 
qu'on execute successivement, a partir du meme jeu de 4p mesures de 
pseudo-distances plusieurs calculs de position du mobile en utilisant 
differentes combinaisons lineaires des phases correspondant aux 

30 frequences L1 et L2. Cela revient k calculer des phases fictives d'une 
frequence porteuse qui serait fictivement la combinaison lineaire des 
frequences L1 et L2. Plusieurs combinaisons lineaires sont successivement 
utilisees pour le meme jeu de 4p pseudo-distances ; le calcul utilise les 
pseudo-distances et une position estinrfee et il aboutit k une position calculee 

35 pour une combinaison Iin6aire donnee ; la position calcufee pour cette 
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combinaison lin£aire sert de position estim^e pour un calcul suivant utilisant 
une autre combinaison lineaire. La succession de combinaisons lineaires, 
correspondant k differentes longueurs d'onde apparentes, est telle que la 
precision du calcul augmente progressivement. Les combinaisons lin£aires 
5 de frequences sont choisies d'un calcul & I'autre notamment de manfere k ce 
que la demi-longueur d'onde apparente correspondante soit supgrieure ou 
6gale k I'erreur de position resultant du calcul precedent, tout en diminuant 
I'erreur ionosph6rique. 

Comme on le verra plus loin, on utilise dans les calculs de 
10 position, des doubles differences 6tablies pour des couples de satellites 
determines de fagon k augmenter encore la longueur d'onde apparente en 
prenant en compte la geometrie de la constellation des satellites (la longueur 
d'onde apparente d'un couple de satellites est d'autant plus grande qu'il est 
vu sous un angle etroit). 
15 En combinant judicieusement le choix des combinaisons lineaires 

de frequences et des couples de satellites, on peut ainsi passer d'une 
position estimee de precision d£cim6trique, voire m6trique, k la position de 
precision centim6trique recherchee. Les combinaisons lineaires sont choisies 
dans une liste dont un exemple est donne plus loin. Cette liste est stockee 
20 par exemple dans I'organe 24 pour etre utilis6e par I'organe de calcul 18. 

Dans le cas du GPS, les combinaisons lineaires s'appliquent aux 
frequences L1 et L2. Le resultat d'une combinaison lineaire de L1 , L2 est 
une nouvelle frequence L3 k laquelle correspond une longueur d'onde dite 
longueur d'onde apparente. Le circuit de reception 16 fournit a I'organe de 
25 calcul 18, la phase <pu du signal 6mis par le satellite pour L1 et la phase <pi_2 
du signal 6mis par le satellite pour L2, ces deux phases 6tant representatives 
de la distance entre le satellite et le mobile ; I'organe de calcul 18 pourra 
alors utiliser non pas la phase cpu ou la phase cpL2 mais une phase cpu* 
correspondant k la combinaison lineaire suivante 
30 L3 = a.L1 + b.L2 selon la formule (pL3 = a. <pu + b. <pl2 

De meme, la station de reference envoie au mobile des phases 
representatives, pour chaque satellite, de la distance entre le satellite et la 
station de reference et la meme combinaison lineaire peut etre utilis6e pour 
determiner une phase apparente k la longueur d'onde L3. 
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Les nombres a et b sont les coefficients de la combinaison lineaire 

choisie. 

Pour chaque combinaison lineaire, on sait calculer non seulement 
la longueur d'onde apparente (done la distance au-dela de laquelle la mesure 
5 de phase devient ambigue) mais aussi la sensibilite aux erreurs de 
propagation ionospherique. Si on prend comme valeur d'erreur 
ionospherique de reference I'erreur existant a la frequence L1, on peut 
calculer la valeur correspondante de I'erreur pour chaque combinaison 
lineaire. Si par exemple I'erreur sur L1 est de 1 cm/km d'eloignement, pour un 
10 eloignement de 10 km a plusieurs dizaines de kilometres, I'erreur sur L1 
atteindrait 10 cm a plusieurs dizaines de cm. 

Pour corriger cette erreur, on pourrait imaginer d'utiliser 
directement une combinaison lineaire de frequences donnant une tres faible 
erreur de propagation ionospherique ; la combinaison dite « lono-Free » 9L1- 
15 7L2 a une tres faible erreur ; mais elle correspond a une longueur d'onde 
apparente tres courte (5cm) de sorte qu'il est impossible de lever son 
ambiguite de phase a partir d'une position estimee eloignee. 

Pour le GPS et a titre d'exemple prefere, on a etabli la liste 
ordonnee suivante de combinaisons lineaires, avec la longueur d'onde 
20 apparente correspondante et le coefficient d'erreur ionospherique 
(relativement a la valeur de reference unitaire pour la frequence L1). Les 
longueurs d'onde apparentes sont decroissantes, et les coefficients d'erreur 
ionospherique decroissent encore plus rapidement. 

L1-L2 connue de I'homme du metier sous le nom de « Wide-Lane », 
25 qui correspond a une longueur d'onde apparente d'environ 86 cm avec un 
ratio de I'erreur ionospherique par rapport a la frequence L1 de base de 1 ,3 ; 

2L1-L2, qui correspond a une longueur d'onde apparente d'environ 16 
cm avec un ratio de I'erreur ionospherique par rapport a la frequence L1 de 
base de 0,56 ; 

30 3L1-2L2, qui correspond a une longueur d'onde apparente d'environ 

13 cm avec un ratio de I'erreur ionospherique par rapport a la frequence L1 
de base de 0,3 ; 

4L2-3L1 , qui correspond a une longueur d'onde apparente d'environ 
1 1 cm avec un ratio de I'erreur ionospherique par rapport a la frequence L1 
35 de base de 0,1 ; 
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9L1-7L2, denommee comme on Pa dej& vu combinaison « lono-Free » 
qui correspond & une longueur d'onde apparente d'environ 5 cm avec un 
ratio de Perreur ionospherique par rapport a la frequence L1 de base proche 
de zero. 

5 Bien sur cette sequence de combinaisons lineaires n'est qu'un 

exemple adapte au GPS dans sa situation actuelle : elie difterera notamment 
dans le cas de frequences L1 et L2 diff6rentes de celles citees (futures 
evolutions du GPS, autres systemes radio-satellitaires tels que Galileo) ; elle 
peut aussi ne contenir que certaines de ces combinaisons lineaires et/ou en 

10 incorporer d'autres. 

Le principe est qu'avec un jeu de 4p pseudo-distances combin6es 
Iin6airement, on peut effectuer un calcul en utilisant une position estimee et 
la premiere combinaison lin§aire de la liste, avec un faible risque d'ambiguTte 
de phase du fait de la longueur d'onde importante de cette combinaison ; on 

15 aboutit a une premiere position calculee plus precise que la position estim6e. 
Cette position calculee sert de position estimee pour un autre calcul fait avec 
les memes 4p pseudo-distances mais combines selon la combinaison 
suivante de la liste, correspondant & une longueur d'onde apparente plus 
faible mais encore suffisamment grande pour ne pas introduire d'ambiguTte 

20 de phase compte-tenu de Perreur ionospherique subsistant suite au premier 
calcul. On aboutit a une nouvelle position calculee toujours plus precise. 
Cette position plus precise sert de position estimee pour un troisfeme calcul 
avec le meme jeu de pseudo-distances combinees selon la troisieme 
combinaison Iin6aire de la liste, et ainsi de suite. A chaque etape la position 

25 est amelior6e du fait de la reduction de Perreur ionospherique, et 
I'amelioration de cette position permet de choisir, dans P6tape suivante et 
sans risque d'erreur d'ambigu!t6, une longueur d'onde plus courte qui avec 
une reduction toujours plus importante de Perreur ionospherique permet de 
rapprocher encore la position calculee de la position vraie. 

30 La figure 3 repr6sente schematiquement les grandes etapes des 

calculs effectu6s. 

Un jeu de 4p pseudo-distances est fourni k Porgane de calcul 
(6tape 1) qui a determine par ailleurs une position estimee Pe et qui choisit 
la premiere combinaison Iin6aire de frequences de la liste predeterminee 

35 (etape 2). 
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Facultativement, le calcul d'une position a partir d'une position 
estimee et d'une combinaison lineaire donnee peut se faire en deux grandes 
etapes : 

- une position approchee Pa est calculee a partir de la position 
5 estimee Pe, de la combinaison lineaire de frequences choisie et d'un jeu de 

4p' pseudo-distances (etape 3a). 

- une position precise Pp est ensuite calculee a partir du point 
approche Pa, de la memo combinaison lineaire de frequences choisie et du 
jeu complet de 4p pseudo-distances (6tape 3b). 

10 Dans une variante de I'invention, on peut calculer directement en 

une etape la position precise Pp a partir de la position estimee Pe, de la 
combinaison lineaire de frequences choisie et d'un jeu complet de 4p 
pseudo-distances ; on peut par exemple retenir le calcul en deux etapes pour 
la premiere combinaison lineaire et en une etape pour les combinaisons 
15 lineaires suivantes. 

Si toutes les combinaisons lineaires de la liste n'ont pas encore 
ete utilisees (§tape 4), on choisit la combinaison lineaire suivante de la liste 
et on reitere les etapes 3a et 3b en considerant comme position estimee la 
position Pp calculee resultant de I'etape 3b (etape 5). 
20 Si toutes les combinaisons lineaires de la liste ont 6t6 utilisees, on 

valide (etape 6) la position Pp resultant du dernier calcul en determinant 
classiquement la valeur de coherence de la position calculee, pour rejeter les 
solutions qui ne respecteraient pas des criteres de coherence minimaux. A 
Tissue de cette etape de validation, la position Pp est la position du mobile 
25 etablie avec une precision centimetrique. 

Dans le cas ou la position Pp n'est pas validee, on execute les 
etapes decrites a partir d'un autre jeu de 4p pseudo-distances. 

On va donner maintenant, a titre d'exemple, le detail du calcul des 
etapes 3a et 3b qu'on peut faire pour Pobtention d'une position precise Pp a 
30 partir d'une position estimee Pe. 

La phase d'obtention d'une position approchee Pa (etape 3a) 
utilise une position estimee initiate Pe ; et elle utilise avantageusement un 
sous-jeu de 4p' pseudo-distances choisies dans le jeu de 4p pseudo- 
distances. 
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On effectue le calcul des positions par un traitement en double 
difference, sur les 4p' pseudo-distances. 

Le traitement dit « en double difference » consiste k travailler non 
pas directement k partir des differences pour deux satellites entre pseudo- 
5 distances mais a partir de differences entre la station de reference et le 
mobile, de la difference pour deux satellites entre pseudo-distances. 

Ces couples de satellites sont choisis en fonction de leur 
sensibility aux erreurs de position, cette sensibility dependant de la 
geometrie entre les satellites et le point de mesure. La g£ometrie de visibility 
10 des p satellites a Pinstant de la mesure etant connue grace aux 6ph6merides 
et k la position estimee, on sait classer les couples de satellites dans Pordre 
de leur sensibility croissante aux erreurs de position. On ne prend que les p' 
satellites correspondant aux couples les moins sensibles, ce qui permet 
d'augmenter la longueur d'onde apparente. Plus precisyment, on calcule les 
15 doubles differences du type 

DDij = (Dim-Djm)-(Dir-Djr), dans lesquelles : 

Dim est la pseudo-distance du mobile au satellite de rang i 
Djm est la pseudo-distance du mobile au satellite de rang j 
Dir est la pseudo-distance de la station de reference au satellite 

20 de rang i 

Djr est la pseudo-distance de la station de reference au satellite 

de rang j 

Les diffyrences du type Dim-Djm ou Dir-Djr permettent d'yiiminer 
les erreurs communes aux satellites ( diffyrences d'horloge entre les 
25 satellites et le r£cepteur). 

Les diffyrences entre ces differences, ou doubles differences DDij, 
permettent d'yiiminer les erreurs dues k la propagation atmospherique ou 
ionosphyrique. 

Les doubles differences sont calculees d'une part pour la 
30 fryquence L1 (difference DDijl), d'autre part pour la fryquence L2 (diffyrence 
DDij2). On calcule alors, pour une premiere combinaison linyaire de 
fryquences, L3 = aL1+bL2, une combinaison lineaire de doubles differences 
CLij = aDDij! +bDDij2. Ce sont ces combinaisons lineaires qui vont etre 
utilisees, et non les classiques doubles diffyrences, pour le calcul de 



WO 03/076962 




PCT/FR03/00746 



position ; elles represented une phase a la longueur cTonde apparente de la 
frequence fictive L3 et sont exprim€es ici en distances. 

Les combinaisons de doubles differences CLij sont comparees a 
des combinaisons similaires, calcutees et non mesurees, k partir de la 
5 position estim6e initiale Pe. La difference qui r6sulte de cette comparaison 
est appelee INNOVij, repr§sentant l'6cart entre r estimation et la mesure. 

Ces ecarts sont relics aux hearts de longitude, latitude, et altitude 
DL, DG et DA entre la position estimee (ici la position initiale Pe) et la 
position calcul6e (ici la position calculee approchee Pa) par des equations 
10 du type 

INNOVij = DL[cos(Evi)cos(Azi)-cos(Evj)cos(Azj)] 
+ DG[cos(Evi)sin(Azi)-cos(Evj)sin(Azj)] 
+DA[sin(Evi)-sin(Evj)] 

ou Evi, Evj sont les elevations des satellites i et j, et Azi, Azj leurs 

15 azimuts. 

Un calcul simple, ou un calcul matriciel avec minimisation des 
erreurs par la technique des moindres carres si on a plus de 4 satellites, 
permet de determiner DL, DG, DA qui represented des hearts entre position 
mesurge et position estimee. Ces ecarts sont ajoutes a la longitude, la 
20 latitude, et Taltitude de la position estimee Pe pour obtenir une position 
approchee Pa. 

A partir de cette position approchee Pa on effectue une deuxfeme 
etape de calcul (etape 3b). La deuxidme etape est tres semblable k la 
premiere, mais 

25 - elle utilise toutes les 4p pseudo-distances, e'est-a-dire celles 

correspondant a tous les couples de satellites, 

- elle utilise comme position estimee non pas la position 
estimee initiale Pe mais la position approchee Pa 

- elle utilise pour le calcul d'une position precise Pp un calcul 
30 matriciel avec un nombre d'6quations en general superieur au nombre 

d'inconnues (le nombre p de satellites £tant suppose superieur a 4) ; la 
determination des hearts DL, DG, DA entre position estimee et position 
calculee peut se faire alors classiquement par une methode des moindres 
carr6s (la position d6termin£e par calcul est cede qui minimise la valeur 
35 quadratique moyenne des residus). 
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La position Pp est ensuite utilis6e dans un nouveau calcul comme 
position estim£e, a la place de la position Pe, et en faisant intervenir la 
combinaison lingaire suivante de la liste. 

Progressivement, par it6ration jusqu'au bout de la liste, on arrive k 
5 une valeur de plus en plus precise et non ambigue de la position du mobile. 
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REVENDICATIONS 

1 . Precede de determination de la position relative d'un mobile par 
rapport a la position connue d'une station de reference, utilisant chacun une 
antenne de reception de signaux radio issus d'une constellation de satellites 
de positionnement 6mettant sur au moins deux frequences L1 et L2, ce 

5 procede comprenant la determination periodique pour chacune desdites 
frequences, d'un jeu de 2p pseudo-distances, a savoir p pseudo-distances 
entre le mobile et les p satellites et p pseudo-distances entre la station de 
reference et les p satellites, la fourniture des pseudo-distances k un organe 
de calcul de position, et le calcul par cet organe d'une position relative du 

10 mobile par rapport a la station de reference k partir d'une part des pseudo- 
distances et d'autre part d'une position estimee Pe du mobile par rapport a la 
station de reference, ce procede etant principalement caracteris6 en ce que 
le calcul de position relative comprend, pour un jeu donne de 4p pseudo- 
distances re?u par I'organe de calcul, les etapes suivantes consistant k : 

15 - a) choisir une combinaison lineaire aL1 +bL2 desdites frequences 

L1 et L2 parmi une liste pred§terminee comportant au moins deux 
combinaisons lineaires de frequences, 

- b) calculer les combinaisons lineaires de pseudo-distances 
correspondant a ladite combinaison lineaire et, a partir de ces combinaisons 

20 lineaires de pseudo-distances et de la position estimee Pe, calculer une 
position relative precise Pp du mobile par rapport a la station de reference, 

- c) choisir dans la liste la combinaison Iin6aire suivante si elle 
existe, et reiterer dans ce cas I'etape b) en considerant comme position 
estimee, ladite position precise Pp, et en utilisant le m§me jeu de 4p pseudo- 

25 distances, pour obtenir une position relative encore plus precise, 

- d) reiterer retape c) pour toutes les combinaisons lineaires de la 

liste. 

2. Procede selon la revendication precedente, caracteris6 en ce 
30 que les combinaisons lineaires de la liste sont determinees de telle sorte 

que, d'un calcul k I'autre, les longueurs d'onde correspondantes diminuent 
progressivement et que la sensibilite aux erreurs ionospheriques diminue 
6galement progressivement. 
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3. ProcedS selon Tune quelconque des revendications 
prec6dentes, caracteris6 en ce que la premiere combinaison lineaire de la 
liste est la combinaison L1-L2 (a=1, b=-1)et/ou la derniere combinaison 
lineaire de la liste est la combinaison 9L1-7L2 (a=9, b=-7), L1 et 12 etant les 

5 frequences d'emission des satellites du systeme GPS. 

4. Procede selon la revendication prec^dente, caracterise en ce 
que les combinaisons lineaires interm6diaires sont de preference les 
suivantes (dans I'ordre) : 

1 0 2L1 -L2 (a=2, b=-1 ); 3L1 -2L2 (a=3, b=-2) ; 4L2-3L1 (a=4, b=-3). 

5. Procede selon Tune quelconque des revendications 
precedentes, caracterise en ce que I'etape b) comprend les deux etapes 
suivantes consistant a : 

15 b1) calculer une position relative approchee Pa du mobile par 

rapport a la station de reference, a partir de la combinaison lineaire choisie, 
de Pe et d'un sous-jeu de 4p' pseudo-distances correspondant a p' satellites 
ou p' est inferieur a p et ou les p' satellites choisis dans la constellation de p 
satellites sont ceux qui ont, compte-tenu de la geometrie actuelle de la 

20 constellation, la plus faible sensibilite a une erreur de position estimee, 

b2) calculer une position relative precise Pp du mobile par rapport 
a la station de reference a partir de ladite combinaison lineaire, de Pa et du 
jeu complet de 4p pseudo-distances. 

25 6. Proc§d6 selon la revendication prec6dente, caracterise en ce 

que les etapes b1) et b2) ne sont realisees que pour la premiere 
combinaison lineaire de la liste, une seule etape faisant intervenir les 4p 
pseudo-distances etant re"alisee pour les autres combinaisons lineaires de la 
liste. 

30 

7. Proc6d6 selon Tune quelconque des revendications 
pr6c$dentes, caracterise en ce que les 2p pseudo-distances entre les 
satellites et la station de reference sont determiners par la station de 
reference et envoy^es par radio au mobile qui comporte alors des moyens 
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de reception pour recevoir ces pseudo-distances et des informations de 
datation de la mesure de ces pseudo-distances. 

8. Dispositif de determination de la position d'un mobile par 

5 rapport a une station de reference, comprenant au moins, dans le mobile, 
des moyens pour recevoir des signaux de positionnement par satellite et des 
moyens pour recevoir un jeu de 2p pseudo-distances emises par la station 
de reference et representant les pseudo-distances entre la station de 
reference et p satellites pour au moins deux frequences porteuses diffe>entes 

10 L1 et 12, des moyens de determination pe>iodique d'un jeu de 2p pseudo- 
distances entre le mobile et les p satellites, des moyens de fourniture des 4p 
pseudo-distances a un organe de calcul (18) de position, des moyens de 
stockage d'une liste de combinaisons lineaires des frequences des porteuses 
des signaux de positionnement, des moyens pour effectuer, a partir d'un 

15 meme jeu de 4p pseudo-distances, des calculs successifs de position 
relative du mobile par rapport a la position de la station de reference, chaque 
fois a partir d'une combinaison lineaire differente de frequences choisie dans 
la liste, d'une position estimee Pe et du jeu de 4p pseudo-distances, la 
position estimee lors d'un calcul avec une combinaison lineaire donnee de la 

20 liste etant la position relative calculee a partir de la combinaison lineaire 
precedente de la liste. 
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